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After decades of painstaking research, the field of heavy ion physics has reached 
, an exciting new era. Evidence is mounting that we can create a high temperature, 

f— *^ ■ high density, strongly interacting "bulk matter" state in the laboratory — perhaps 

^/■^ ' even a quark-gluon plasma. This strongly interacting matter is likely to provide 

f— *^ ' qualitative new information about the fundamental strong interaction, described 

""-"^ , by Quantum Chromodynamics (QCD). These lectures provide a summary of ex- 

^ i perimental heavy ion research, with particular emphasis on recent results from 

RHIC (Rclativistic Heavy Ion Collider) at Brookhaven National Laboratory. In 
addition, we will discuss what has been learned so far and the outstanding puzzles. 

O ' 
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1. Introduction 

While the universe as we know it is well described by the standard model of 
particle physics, some important questions remain unanswered. Perturba- 
tive Quantum Chromodynamics (pQCD) — a part of the standard model 
— is a very successful description of hard, or short-distance, phenomena , 
where the "strong interaction" becomes weak due to asymptotic freedom. 
For example, the production of jets in pp collisions at 1.8 TeV is well de- 
scribed for jet transverse energies from 10-400 GeV 2 . There is, however, 
an important set of soft physics phenomena that are not well understood 
from first principles in QCD: color confinement, chiral symmetry breaking, 
and the structure of the vacuum. These phenomena are important: almost 
all of the visible mass of the universe is generated by soft QCD and not 
by the direct Higgs mechanism. The current masses of the three valence 
quarks make up only about 1% of the mass of the nucleon 3 . 

In order to study these phenomenon, we seek to separate color charges 
by heating matter until a quark-gluon plasma is formed. A conventional 
electromagnetic plasma occurs at temperatures of about 10 4 -10 5 K, cor- 
responding to the typical ionization energy scale of 1-10 eV. Theoretical 
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studies of QCD on the lattice 4 indicate that the typical energy scales of 
thermally driven color deconfincment are in the vicinity of 170 MeV, or 
2 x 10 12 K. In addition to providing information about the strong inter- 
action, achieving such temperatures would also provide a window back in 
time. The color confinement phase transition is believed to have occurred 
within the first few microseconds after the big bang. 

In order to achieve such high temperatures under laboratory conditions, 
it is necessary to use a small, dynamic system. For instance, experimental 
fusion reactors heat a conventional plasma up to temperatures as high as 
10 8 K over distance scales of meters and lasting for seconds. By colliding 
gold ions at nearly the speed of light, we expect to achieve temperatures of 
order 10 12 K over distance scales of order 10 fm and time scales of order 10- 
100 ys a . Clearly one of the challenges in this endeavor will be to determine 
whether such small and rapidly evolving systems can elucidate the bulk 
behavior that we are interested in. Another challenge will be to use some 
of the rarer products of the collisions to probe the created "bulk" medium. 

The focus of these lectures will be on the results coming out of the Rel- 
ativistic Heavy Ion Collider (RHIC) experiments at Brookhaven National 
Laboratory (BNL). Earlier experimental results and some theoretical work 
will be mentioned as needed, but a comprehensive review of heavy ion 
physics will not be attempted. The RHIC spin physics program using po- 
larized protons will also not be covered. 

2. The Machine and Detectors 

The RHIC data described in these lectures were taken during the last three 
years of running at RHIC, starting in the summer of 2000, as summarized 
in Table 1. The runs were characterized by their species and their y/s NN , 
which is the cm collision energy of one nucleon taken from each nucleus. 
For instance, a AuAu collision with 100 x A GeV on 100 x A GeV would 
have y/s NN = 200 GeV. Most of the runs were several weeks in duration, 
with two exceptions. The 56 GeV run, not intended as a physics run, was 
only 3 hours long and data is only available from a preliminary subsystem 
of one experiment (PHOBOS). The 19.6 GeV run was 24 hours long and 
usable data were taken by three experiments. For the 130 and 200 GeV 
runs, all four detectors participated: two large detectors/collaborations 
with 300-400 collaborators each and two small detectors/collaborations 



a Rccall that one yoctosecond = 10 24 s. 
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with 50-70 collaborators each. These four detectors complement each other 
and have provided a broad range of physics results. The BRAHMS exper- 
iment (Broad RAnge Hadron Magnetic Spectrometer) focuses on tracking 
and particle ID at high transverse momentum over a broad range of rapid- 
ity from 0-3. The PHENIX experiment (Pioneering High Energy Nuclear 
Interaction experiment) provides a window primarily at mid-rapidity, but 
specializes in high rate and sophisticated triggering along with a capability 
to measure leptons and photons as well as hadrons. The PHOBOS ex- 
periment (descendant of the earlier MARS experiment) provides nearly Att 
coverage for charged particle detection, good vertex resolution, and sensi- 
tivity to very low p T particles. The STAR experiment (Solcnoidal Tracker 
At Rhic) provides large solid angle tracking and complete coverage of ev- 
ery event written to tape. More details concerning the capabilities of the 
accelerator and experiments can be found in NIM journal issue dedicated 
to the RHIC accelerator and detectors 5 . 



Table 1. RHIC running conditions to date. 



Date 


Species v^T^" 


Year 2000 
Year 2001 
Jan. 2002 
Year 2003 


AuAu 56, 130 GeV 
AuAu 19.6, 200 GeV 
pp 200 GeV 
dAu, pp 200 GeV 



Some data will also be shown from lower energy heavy ion collisions, 
particularly from the CERN-SPS (Conscil European pour la Recherche Nu- 
clearie - Super Proton Synchrotron) will also be discussed. The top CERN 
energy is ^/s^ =17.2 GeV. 

3. Strongly Interacting Bulk Matter 

In order to learn anything about QCD from heavy ion collisions, we must 
first establish that we have created a state of strongly interacting bulk 
matter under extreme conditions of temperature and pressure. 

3.1. How Much Matter? 

Figure 1 shows the charged particle distribution for central (head-on) AuAu 
collisions in the pscudorapidity variable: r\ = — lntan(#/2). These data 
imply a total charged multiplicity of 1680 ± 100 for the 19.6 GeV data and 
5060 ± 250 for the 200 GeV data 6 . While this number is considerably 
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Figure 1. Pseudorapidity distributions, dN^/dr), for central (6%) AuAu collisions at 
V /577^T = 19.6, 130, and 200 GeV. Data taken from PHOBOS 6 . 

smaller than Avagadro's number, it is substantial thermodynamically since 
small-system corrections to conventional thermodynamics start to become 
unimportant for systems with about 1000 particles or more 7 . 

The number of particles produced in a given AuAu collision varies widely 
due to the variable geometry of the collision. Some collisions are nearly 
head-on with a small impact parameter, while most collisions have a larger 
impact parameter, with only a partial overlap of the nuclei. These cases 
can be sorted out experimentally, using both the number of produced par- 
ticles and the amount of "spectator" neutrons seen at nearly zero degrees 
along the beam axis. The impact parameter or "centrality" of the collision 
is characterized by the number of nucleons from the original ions which 
participate in the heavy ion collision, (N part ), or the number of binary NN 
collisions, (N co u). More details can be found in Refs. 8 ' 9 . 

3.2. Elliptic Flow: Evidence for Collective Motion 

Non-central heavy ion collisions have an inherent azimuthal asymmetry. 
The overlap region of two nuclei is roughly ellipsoidal in shape. If there is 
collective motion that develops early in the collision, this spatial anisotropy 
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Figure 2. Left panel: elliptic flow as a function of centrality as seen by STAR (data) 
compared to hydrodynamic models (rectangles) 10 . Right panel: peak elliptic flow as a 
function of collision energy for ultrarelativistic collisions, taken from an NA49 compila- 
tion 11 . 



can be converted to an azimuthal asymmetry in the momentum of detected 
particles. This azimuthal asymmetry is characterized by a Fourier decom- 
position of the azimuthal distribution: 

dN/d<p = N (l + 2vi cos (f> + 2v 2 cos(20)), (1) 

where cj> is the azimuthal angle with respect to the reaction plane b . The 
left-hand panel of Fig. 2 shows that the elliptic flow parameter is quite 
large, nearly reaching the values predicted by hydrodynamic models. These 
models assume a limit of local equilibrium with collective motion of the bulk 
"fluid" . The right-hand panel of Fig. 2 shows that this asymmetry is the 
largest ever seen at relativistic energies. 

Elliptic flow, in addition to indicating that there is collective motion, 
can provide information about the type of motion. In particular, the p T de- 
pendence of elliptic flow can distinguish between two limits: the low density 
limit and the hydrodynamic limit (rapidly expanding opaque source). In 
the low density limit, some of the produced particles are absorbed or scat- 
tered once (and usually only once). In this case, for relativistic particles, 
V2 is nearly independent of p T . In the hydrodynamic limit, in contrast, we 
expect t>2 oc p T for moderate values of p T . This effect comes about because 
the expansion causes a correlation between normal space and momentum 



b The true reaction plane is defined by the impact parameter vector between the gold ions. 
The experimental results shown have been corrected for the reaction-plane resolution, 
which would otherwise dilute the signal. 
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Figure 3. Elliptic flow versus p T for all particles (left panel) 10 , and for identified parti- 
cles (right panel) from STAR 12 . The curves in the right panel refer to a hydrodynamic 
model description. 



space, forcing the highest p T particles to come from the surface, while low 
p T particles can come from anywhere in the volume. Data from the SPS 
favor the hydrodynamic limit 13 . The left-hand panel of Fig. 3 shows a 
clear linear relationship between elliptic flow and transverse momentum at 
RHIC as well, while the right-hand panel shows that hydrodynamic models 
not only describe the overall trend, but even describe the pions and protons 
separately. 
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Figure 4. Elliptic flow as a function of pseudorapidity from PHOBOS 14 . 



Finally, elliptic flow can be examined as a function of pseudorapidity. 
The expectation was that the elliptic flow would be nearly independent of 
pseudorapidity as the basic physics of RHIC were expected to be invariant 
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under longitudinal boosts. Fig. 4 shows that t>2 is strongly dependent on 
pscudorapidity, a result which has still not been explained. 

Taken together, these results show clear evidence of collective motion 
and suggest a system at or near hydrodynamic equilibrium which is rapidly 
expanding in the transverse direction and which does not exhibit longitu- 
dinal boost-invariancc. 

3.3. Hanbury-Brown Twiss Effect: More Dynamics 

Intensity intcrfcrometery, or the Hanbury-Brown Twiss effect 15 , is a tech- 
nique used to measure the size of an object which is emitting bosons (e.g. 
photons from a star or pions from a heavy ion collision). Boson pairs which 
are close in both momentum and position are quantum mechanically en- 
hanced relative to uncorrelated boson pairs. Bosons emitted from a smaller 
spatial source are correlated over a broader range in relative momentum, 
which allows you to image a static source using momentum correlations. 

For a given pair of identical particles, we can define their momentum 
difference, q, and their momentum average, k. We can further define the 
three directions of our coordinate system 16 : 

• Longitudinal (Ri) — along the beam direction (i), 

• Outwards (R D ) — In the (z, k) plane, _L z, 

• Sidewards (R s ) — _L z & _L k. 

For a boost-invariant source, the measured sidewards radius at low p T 
will correspond to the actual physical transverse (rms) extent of the source 
at freezeout, while the outwards radius will contain a mixture of the spatial 
and time extent of the source. Particles emitted earlier look like they are 
closer to the observer, which artificially extends the apparent source in the 
out direction. In particular, 



of emission" parameter, a x and a y are the geometric size in the out and 
side directions, and a XT is the space-time correlation in the out direction. 

In the case of an azimuthally symmetric and transparent source, the 
last two terms are taken to be small or zero and we have 



Given the assumption of a boost-invariant, azimuthally symmetric and 
transparent source, the HBT results from heavy ion collisions have been 



where 0]_ 




R — R 2 S ~ 0±.°t- 



(3) 
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Figure 5. HBT parameters as a function of colliding beam energy. 



perennially confusing. From Eq. 3, we expect R D /R S > y2 since most 
sources should emit for a time which is of the same order as their size. 
Some models of the Quark-Gluon Plasma predict an even larger value for 
this ratio as the plasma might need to emit particles over a long time du- 
ration in order to get rid of the entropy 18 . However, as can be seen in 
Fig. 5, R /R s is basically unity at RHIC energies, naively implying an 
instantaneous emission of particles over a moderately large volume. 

This situation, along with the modest values of Ri , has been termed the 
"HBT puzzle" . Primarily, though, these data indicate a need to improve the 
modeling of the collision. If you consider a source which is opaque, rapidly 
expanding and also not boost invariant, the meaning of — R 2 S changes 
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since we must use Eq. 2 and not Eq. 3. Opacity reduces the apparent R D 
value since you only see the part of the source closest to you in the out 
direction. Transverse expansion along with opacity will decrease the ratio 
further since particles emitted later are also emitted closer to the viewer, 
reducing the magnitude of R . Finally, a general longitudinal expansion 
(not just coasting) must be taken into account since we know that the 
source is not boost invariant. This effect would explain the small size of Ri 
and has also been shown 19 to reduce the ratio R /R a . 

So, while HBT and elliptic flow have not been successfully described 
in full detail by the hydrodynamic models yet, the qualitative message 
they provide is very similar. The source is rapidly expanding (probably in 
all three dimensions), opaque, and can be described as "hydrodynamically 
equilibrated bulk matter" . 

3.4. Characterizing the Bulk Matter 

Having established that the system has a large number of particles as well 
as collective behavior, we can now proceed to consider bulk quantities such 
as the temperature and baryon chemical potential of the system. 

In conventional, static, thermodynamic systems, the temperature can 
be measured by directly measuring the average energy per particle. In a 
very dynamic system, such as a heavy ion collision, we have to separate 
the energy contributed by collective motion from the thermal energy. To 
do this, we make use of the fact that the collective velocity contributes 
more to the momentum of heavy particles than to lighter particles. Ther- 
mal fits 20 to (p T ){m) yield a temperature of approximately 100 MeV and 
an average transverse expansion velocity of 0.55 c. This large expansion 
velocity supports the picture given by the elliptic flow and HBT. 

Another thermometer is provided by the fact that the ratios of particles 
of different masses are sensitive to the temperature. In addition, ratios 
of particles with the same mass, but different quark content, such as p/p 
and K~ /K + , are sensitive to the balance between matter and antimatter, 
characterized by the baryon chemical potential fis- Positive values of fis 
refer to a matter (baryon) excess in a system. Fig. 6 shows particle abun- 
dance ratios and a thermal fit. This fit yields a constant temperature of 
176-177 MeV at both energies, but a falling baryon density (41 MeV at 
130 GeV and 29 MeV at 200 GeV). The falling baryon density is expected. 
Higher energy collisions dilute the fixed initial baryon excess from the orig- 
inal gold nuclei and also make it harder to transport the baryon excess to 
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Figure 6. Particle abundance ratios and a thermal fit for the two highest RHIC ener- 
gies 21 . 



midrapidity. 

We are immediately faced with a dilemma: our kinetic thermometer, 
based on energy per particle, indicated a temperature of ~ 100 MeV, while 
our chemical thermometer, based on particle abundances, indicated a much 
higher temperature of ~ 175 MeV. The resolution of this paradox lies in the 
fact that only inelastic collisions can change the particle abundances while 
both elastic and inelastic collision serve to equalize the energy between par- 
ticles. Using the terminology of cosmology, we can define an approximate 
"freezeout hypersurface" which contains the spacetime points of the final 
collisions suffered by each particle. In the case of a heavy ion collision, 
the chemical freezeout can occur earlier than the kinetic freezeout. This 
resolves our dilemma, but with the unavoidable consequence of making our 
picture of the collision somewhat more complicated. It should be noted 
that the HBT results are actually imaging the kinetic freezeout boundary 
as the source. 



3.5. "Little Bang Cosmology" 

As in cosmology, we are interested in understanding what happened before 
the freezeout. We can estimate the energy density from the transverse 
energy produced in the collision and the cylindrical volume occupied shortly 
after the collision occurred. This leads to the formula 22 : 

1 1 dE T 

e = (4) 

7ri? 2 ctq dy 

where the radius R is the nuclear radius and To is the time it takes for the 
transverse energy to be effectively equilibrated (0.2-1.0 fm/c). 
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Figure 7. Transverse energy at midrapidity as a function of centrality for 130 GcV 
and 17 GeV collisions. Left panel: per participating nucleon, right panel: per produced 
particle. Data taken from 23 . 




Figure 8. Charged particle multiplicity per participating nucleon pair at midrapidity 
as a function of beam energy 24 . 

Figure 7 (right panel) shows that the transverse energy per particle is 
about 800 MeV at RHIC while Fig. 8 shows the multiplicity. Combining 
these results, using Eq. 4 yields e =5-25 GeV/fm for central collisions at 
the highest RHIC energy. Fig. 9 shows the theoretical relationship, based 
on lattice QCD calculations 4 , between energy density and temperature. 
The expected T 4 dependence of an ideal gas has been divided out, leading 
to a a constant value for high temperature, proportional to the number of 
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Figure 9. Energy density scaled by T 4 (natural units h = c = k = 1) as as a function of 
temperature scaled by the critical temperature (T/T c ). The arrows on the right indicate 
the Stefan-Boltzmann values for an ideal non-interacting gas. Figure taken from Ref. 4 . 

degrees of freedom in the quark-gluon plasma phase. 

Combining the data with the theoretical curves leads to an estimated 
initial temperature of 300 ± 50 MeV for central AuAu collisions at the top 
RHIC energy. This is significantly higher than the theoretical transition 
temperature of ~ 170 MeV. A similar exercise at the top CERN-SPS energy 
y/s NN =17 GeV, yields an estimated initial temperature of 240 ± 50 MeV. 
It should be noted that if we assume a hadronic description rather than 
a phase transition, the number of degrees of freedom should actually be 
lower, implying an even higher initial temperature (about twice as high). 
This means that the estimated initial temperatures of ~ 300 and 240 MeV 
for RHIC and CERN actually represent lower limits. 

3.6. Summary: Bulk Matter 

Figure 10 shows the phase diagram based on the chemical freezeout points 
measured at various energies in heavy ion collisions including the 130 GeV 
point from RHIC. The 200 GeV point from RHIC would be at basically the 
same temperature, but \xb = 29 MeV rather than fis = 41 MeV. The curve 
through the data implies freezeout at a fixed energy per particle of about 
1 GeV, while the bands indicate the theoretical expectation for the tran- 
sition between confined and deconfined matter. The initial temperatures 
estimated for both RHIC and CERN are not shown, but they would lie 
above the theory curve, with the RHIC temperature being 300 ± 50 MeV. 
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Figure 10. Phase diagram of heavy ion collisions from Ref. 25 . The data points represent 
heavy ion collisions over a broad range of energies. The curve through the data points 
represents a fixed energy per particle. The upper band represents an estimate of the 
phase boundary. The lower band represents a constant energy density (0.6 GeV/fm 3 ). 
The isolated point above the theory curves represents a theoretical critical point. 



The constant freezeout temperature for high energy ion collisions, appear- 
ing at the theoretical boundary between confined and deconfined matter is 
provocative. It could be an accident, but it is similar to a situation where 
you have a detector which only detects liquid, you determine indirectly that 
you created matter at 200° C, and you directly detect droplets of water at 
a temperature of 100°C. 

To summarize this section, we have produced a dense, hot, rapidly ex- 
panding bulk matter state. We have seen a universal freezeout curve and 
it is suggestively close to the expected boundary between deconfined and 
confined matter. Furthermore, we have indications that the initial collision 
reaches energy densities (and therefore temperatures) well in excess of that 
expected to be needed for deconfincmcnt. 

Efforts to probe this state quantitatively are just beginning, but show 
promise. This will be the subject of the next section. 



February 8, 2008 21:17 WSPC/Trim Size: 9in x 6in for Proceedings 



baker'LL03 



14 

4. Probing the Earliest, Hottest Part of the Collision 

While the freezeout temperature measurements are on solid footing, the 
estimates presented above for the initial temperature are indirect. As the 
RHIC program develops, we can go beyond these qualitative discussions of 
the early times and start probing them more quantitatively. 



4.1. Electromagnetic and Hidden Flavor Probes 

Perhaps the cleanest method, from a theoretical perspective, would be to 
examine thermal photons and leptons that originate from the early part of 
the collision when the temperature was higher. These weakly interacting 
particles are expected to decouple thermally (or "freezeout" ) from the bulk 
strongly interacting matter much earlier than hadrons. Combined with 
a measurement of the energy density this would effectively measure the 
number of degrees of freedom in the initial state. While it is theoretically 
very clean, this measurement is experimentally very challenging. A typical 
central collision at RHIC produces thousands of neutral pions which decay 
into thousands of photons in each event and serve as a background to this 
measurement. A typical RHIC detector also has literally tons of material 
in which background photons (and leptons) can be produced. 

Despite the difficulty, these measurements and fits have been attempted 
at the SPS, both in terms of direct photon spectra 26 and thermal leptons 27 . 
These results lead to an estimated initial temperature at the SPS of ~ 
200 MeV, consistent with our estimate above for partonic matter. These 
results, however, are very sensitive to details of how the backgrounds are 
handled. 

Hidden heavy flavor measurements (strangeness, charm, and beauty) 
also show promise as potentially sensitive probes of the density of the 
medium and of chiral symmetry restoration. Fig. 11 shows the suppres- 
sion of the J/V> (hidden charm) compared to collision scaling at RHIC and 
the SPS. Sensitivity can also be found in the mass, line shape, and yield 
of the 4> particle (hidden strangeness), seen by its hadronic and leptonic 
decay modes 28 . So far at RHIC, these measurements suffer from lack of 
statistical power. 

One common denominator that would make many of these signals 
clearer would be a clean measurement of open heavy flavor (D and B par- 
ticles). These measurements should be forthcoming from RHIC following 
upgrades to the detectors and improved luminosity from the collider. 
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Figure 11. Measured yield of J/ip by NA50 at CERN (stars) and by PHENIX at RHIC 
(black squares and arrows) 29 . 



4.2. Hadron Suppression: Jet Quenching? 

In addition to measuring the initial temperature, we would like to have 
a more direct measure of the energy density of the bulk matter that we 
have created. One handle on this quantity is to study the behavior of high 
momentum particles in heavy ion collisions. In particular, partons with 
relatively high transverse momentum are predicted to lose energy when 
traveling through dense matter, in a phenomenon known as "jet quenching" . 
The amount of energy loss is proportional to the energy density of the 
matter traversed, so this is potentially a very sensitive probe. 

All four experiments at RHIC measured particle spectra 30 > 31 > 32 i 33 > 34 > 35 . 
These spectra need to be compared to a reference sample, appropriately 
scaled. The simplest such reference sample is to consider each NN collision 
in the initial A A collision geometry as being independent. This leads us to 
define a "nuclear modification factor" : 



Raa{pi 



AA 



/dp T dr] 



N, 



coll 



d 2 aPP/dp T dr) ' 



(5) 



At high momentum, (px > 2 GeV/c) this ratio should approach unity if 
the collisions are independent and the jets are not affected by the material. 
Jet energy loss in the medium should show up as a suppression of high 
momentum hadrons. In lower energy AA collisions 36 and pA collisions 37 , 
an excess has been observed rather than a suppression. This effect is inter- 
preted as being caused by multiple scattering during the initial collision. 
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Figure 12. Yields from PHENIX for a) charged particles and b) neutral pions for central 
and peripheral 130 GcV AuAu collisions compared to scaled reference samples 35 . 

The results at RHIC energies are strikingly different from lower energy 
data as can be seen in Fig. 12 from PHENIX. Invariant yields for produced 
particles in central and peripheral 130 GeV AuAu data are compared to a 
scaled-up pp reference sample. For p T > 2 GeV, the peripheral data scales 
as expected, while the central AuAu data shows a substantial suppression. 
The dramatic difference between the different energies is even more appar- 
ent in Fig. 13 where the scaled reference data are divided out to yield Raa, 
the nuclear modification factor of Eq. 5. Clearly, something qualitatively 
different is occurring at RHIC energies. Similar results were seen by STAR 
at 130 GeV 30 and all four experiments at 200 GeV 3i,32,33,34_ Since thc 
peripheral data scales as expected, it is also possible to measure hadron 
suppression by taking the ratio of central/peripheral data, scaled by the 
ratio of N co u. At 200 GcV, this technique was used to establish that this 
hadron suppression persists to very high transverse momentum, as seen in 
Fig. 14. As indicated above, this hadron suppression may be a signature of 
jet quenching, in which case we have clear evidence of a system with very 
high energy density. 

Another view of this hadron suppression, from PHOBOS and PHENIX, 
shows how strong thc effect is. Fig. 15 (left panel) shows the yield in AuAu 
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Figure 13. Nuclear modification factor for charged particles and neutral pions for 
17 GeV PbPb collisions at CERN and central 130 GeV AuAu collisions from PHENIX 35 . 
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Figure 14. Scaled ratio of charged particle production in central to peripheral AuAu 
collisions at 200 GeV from STAR, using two different choices of peripheral data 33 . 
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Figure 15. Scaling of charged hadron yields from AuAu collisions from PHOBOS and 
PHENIX. Left panel: Yield scaled per participant pair ((N par t) /2) normalized to the 
yield of the most peripheral bin (45—50%) 31 . Right panel (lower plot): Yield scaled per 
participant pair with reference to pp data 38 . 



collisions for fixed values of p T , scaled by mid-central data ((N part ) ~ 
65) and normalized by (N part ). Consider the lower-right hand plot, with 
p T — 4.25 GeV/c. The dashed curve shows the expectation if N co U (A 4 / 3 ) 
scaling held true over the centrality range shown, while the solid line shows 
the expectation of N part (A 1 ) scaling . For N part > 65, we see approximate 
A 1 scaling of high p T particle production. Fig. 15 (right panel) shows a 
similar result, normalized to pp data, from PHENIX. The lower plot is 
normalized per participant and the result for p T > 3.6 GeV/c is relatively 
flat for (N part ) > 80. This particular form of high p T suppression could be 
an indication that jet quenching reaches a geometric maximum involving 
one power of length scale Ra u N^ t (see e.g. Ref. 39 ). 

Another piece of evidence in favor of the jet quenching interpretation for 
this data comes from STAR. Jets in pp collisions can be seen by triggering 
on a high momentum particle and then looking for correlations of moderate 
p T particles azimuthally. In pp collisions, this leads to a clear two-jet signal 
with a cluster of particles near the trigger particle in azimuth and another 
cluster at = it (back-to-back correlation). This signal indicates that 
jets are created and acquire large transverse momentum in conventional 
2 — > 2 parton scattering processes and that the jets survive. For peripheral 
AuAu collisions, one expects a similar result as found in pp, with a small 
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correction due to correlations induced by elliptic flow . 



p+p data + flow 



6 < p T (trig) < 8 GeV/c 
• 0-5% Au+Au 

flow: v2 = 7.4% 

pp data + flow 




A o (radians) 

Figure 16. Azimuthal correlation functions from STAR 40 . Left panel: Data from 
peripheral AuAu collisions (filled circles) compared to a jet+flow reference sample (open 
circles) and a flow-only reference curve. Right panel: Data from central AuAu collisions 
(points) compared to a jet+flow reference sample (upper curve) and a flow-only reference 
curve (lower curve). 



Figure 16 (left panel) shows that the reference sample constructed from 
pp collisions and the measured elliptic flow successfully describes the pe- 
ripheral AuAu data: jets are created back-to-back and survive. In contrast, 
Fig. 16 (right panel) shows the result for central AuAu data. In this case, 
the azimuthal correlation function agrees with the jet+flow reference for 
Acj) ~ while it agrees with the flow-only reference for A<fi ~ n. This 
means that the near-side jet survives, but the away-side jet disappears. 
The main point here is that this measurement shows that the hadron sup- 
pression is a jet phenomenon. If back-to-back jets are indeed produced as 
expected in central AuAu collisions, then the away-side jet is quenched by 
the bulk matter. 



4.3. Is Jet Quenching the Only Possible Explanation? 

Triggered by the observation that the scaling is approximately proportional 
to N part , or A , at large p T , Kharzeev, Levin, and McLerran showed that 
the "suppression" of jets in AuAu compared to NN could simply be due to 



c Since particles are preferentially produced in the event plane, a trigger particle in the 
event plane will tend to pick up particles at A</> = or 7r. This means that the appropriate 
reference is C2(p + p) + A(l + 1v\ cos(2A0)). 
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initial state effects already present in the gold nuclear wavefunction 41 . Par- 
ton recombination (or saturation) can cause gluons from different nucleons 
in the gold nuclei to recombine, leading to a smaller number of partons with 
a higher transverse momentum per parton. Qualitatively, this is difficult 
to distinguish from jet quenching since it reproduces both effects: 

(f ) There are fewer high p T jets than expected because the gold nuclei 
are not simple linear superpositions of nucleons and there are just 
fewer quarks and gluons to begin with than expected. 

(2) Jets do not necessarily come out back-to-back. The usual argument 
for back-to-back jets assumes two incoming partons with p T ~ 
followed by a large angle 2 — > 2 scatter into two back-to-back jets. 
However, multiple parton collisions in the initial state lead to par- 
tons with non-zero p T compensated by multiple partners, which 
need not appear at Acf> = ir. 




These authors also showed that the initial state saturation model could 
be made to agree quantitatively with the data (see Fig. 17), including the 
effect of approximate N part scaling 41 . 
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The saturation model also describes the overall particle production. In 
fact, this ability is more natural since parton saturation effects are strongest 
for the softest partons where the parton densities are the highest. The 
saturation model relates the gluon distribution at low x in deep inelas- 
tic scattering from protons with the energy, centrality, and pseudorapidity 
dependence of particle production in heavy ion collisions 42 . 




Figure 18. The charged particle pseudorapidity distributions and Kharzeev-Levin sat- 
uration model fits. Left panel: PHOBOS data at 130 GeV 42 . Right panel: BRAHMS 
data at 200 GeV 43 . 

It should also be noted that the saturation model was one of the few 
models to correctly predict all of the following: the 130 and 200 GeV midra- 
pidity multiplicity 24 > 44 and the centrality dependence at all three ener- 
gies 45 . Figure 18 shows the fits to 130 and 200 GeV data from PHOBOS 
and BRAHMS respectively. The pseudorapidity and energy dependence 
are primarily controlled by the A parameter, which is extracted from deep 
inelastic scattering data. 

So the initial state saturation model describes well the bulk of soft 
particle production and, if pushed, may also describe the moderately high 
p T particle production behavior. More importantly, the hadron suppression 
or "jet quenching" effect which we want to use as a probe of the density of 
the strongly interacting bulk medium may not be a final state effect at all, 
but may be actually be present in the gold wavefunction. 

4.4. Initial or Final State Effect? 

At the time of these lectures, RHIC was running deuteron-gold collisions in 
order to resolve this issue. Initial state effects, such as parton saturation, 
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should still occur in dAu since they are associated with the gold nucleus 
itself and not the collision. Final state effects, such as jet quenching, should 
go away in dAu since we do not expect a large bulk medium to form. Some 
preliminary hints already indicated that the suppression was probably a 
final state effect rather than an initial state effect. 

Since charm quarks are primarily formed by gluon-gluon fusion and are 
not expected to be quenched in the final state 46 , charm serves as a measure 
of the number of gluons available for hard scattering from the initial state. 
Open charm production, which was found to scale with the number of 
collisions 47 , implies that parton saturation does not affect hard scattering. 
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Figure 19. Elliptic flow of charged particles in 130 and 200 GeV AuAu collisions (and 
their ratio) from STAR 48 . 



The behavior of elliptic flow at high p T also suggests that high p T par- 
ticles are strongly absorbed in the final state. Figure 19 shows that elliptic 
flow reaches a constant value at high p T , independent of y/s NN . Further- 
more, the value is so large that it is essentially the maximum allowable 
asymmetry from a geometric point of view 49 . This implies that only jets 
emitted close to the surface make it out as was also indicated by the ap- 
proximate N part scaling of high p T particles. Since the transverse geometry 
of the collision is a final state effect and not present in a single initial gold 
wavefunction, high p T particles must be strongly absorbed or rescattered 
in the final state, such that the collision geometry leaves its imprint on the 
final state momentum distribution. It should be noted that a t>2 value of 
0.17 implies that twice as many particles are emitted in-plane as out-of- 
planc, a huge effect. This effect has been shown by STAR 48 to persist to 
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p T > 8 GeV/c. 

Taken together, the results on overall particle production and high p T 
particle production are very suggestive. The system appears to be made 
up of hydrodynamic bulk matter. The system is opaque and expanding ex- 
plosively, probably in all three dimensions. The estimated energy density 
is much higher than that of the theoretical transition. There is a freezeout 
along a universal curve near the theoretical transition. There is a strong 
suppression of inclusive high p T yields and back-to-back pairs and an az- 
imuthal anisotropy at high p T . The natural implication is that there is a 
large parton energy loss and surface emission. 

These results are tantalizing, but there are some caveats. First of all, we 
not yet have a complete 3D hydrodynamic description of the collision which 
is consistent with all of the data. Additionally, there are some outstanding 
puzzles from PHOBOS and PHENIX. Finally, data from dAu collisions are 
needed to really disentangle initial state effects. We will turn to the puzzles 
and dAu data next. 

5. Some Puzzles at RHIC 

In addition to the surprising features mentioned above (blackness and 3D 
explosiveness of the source), there are two deep puzzles in the data: the 
behavior of protons at moderately high p T and the apparent universality 
of particle production at high ^Js NN . 

5.1. Scaling Puzzle I: Bary on/ Meson Differences 




p T [GeV/c] p T [GeV/c| 

Figure 20. Proton to pion ratio for 200 GeV AuAu collisions from PHENIX 50 . 

The first puzzle, emphasized initially by PHENIX, concerns the remark- 
able number of protons (compared to pions) at large transverse momentum, 
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as shown in Fig. 20. PHENIX has also shown that pions are more sup- 
pressed than protons in the intermediate p T region from 2-5 GeV/c 51 . 
This effect is also seen in neutral mesons and baryons by STAR 52 . 

Why are pions more suppressed than protons? The current ideas in- 
clude a modification to the fragmentation function in the hot medium or a 
difference in gluon jet and quark jet quenching in the hot medium. Perhaps 
the most intriguing explanation is that, in the presence of jet quenching, a 
different production mechanism — quark coalescence — starts to dominate 
hadron production. Instead of the usual jet fragmentation, this is a multi- 
parton mechanism: three independent quarks coalesce into a baryon or an 
independent quark and antiquark coalesce into a meson 53 . 
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Figure 21. Elliptic flow per constituent quark as a function of transverse momentum 
per constituent quark for lambdas and neutral kaons for 200 GeV AuAu 52 . 



Since (1 + 2v 2 cos20) Ar w (1 + 2Av 2 cos20), the coalescence model 53 
predicts a scaling in elliptic flow per constituent quark versus p T per con- 
stituent quark. Figure 21, from STAR, shows this scaling effect for elliptic 
flow. This model also explains the fact that, in AuAu collisions at high 
p T , baryons and mesons behave similarly while mesons are suppressed (and 
reach maximum V2) at lower momentum. 

While this explanation is intriguing, this result remains a puzzle because 
it is unclear that this model should apply to dAu data (see Section 6). 
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Figure 22. Comparison of the total charged multiplicity versus collision energies for 
AA, c+e — , pp, and pp data, as described in the text, from PHOBOS 54 . In the upper 
panel, the curve is a perturbative QCD expression fit to the e+e~ data. In the lower 
panel, the data have all been divided by the e+c~ fit. 



5.2. Scaling Puzzle II: Similarity of A A and e+e at High 
Energy 

Figure 22 shows the total charged multiplicity for AA collisions (scaled by 
(A r par t)/2) compared to pp, pp, and e + e~, as a function of the appropriate 
yfs for each system 54 . The e + e~ data serve as a reference, describing the 
behavior of a simple color dipole system with a large y/s. The curve is a de- 
scription of the c + c~ data, given by the functional form: Ca s (s) /1 eV ,B / a »( s ' 
with the parameters A and B calculable in perturbative QCD and the con- 
stant parameter C determined by a fit to the e + e _ data 55 . In order to 
compare them with e + e~, the pp and pp data were plotted at an effec- 
tive energy ^/s c ff = which accounts for the leading particle effect 56 . 
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Finally, central AA collisions, AuAu from the AGS and RHIC, and PbPb 
from CERN are shown. Over the available range of RHIC energies from 
19.6 to 200 GeV, the AuAu results are consistent with the e + e~ results, 
suggesting a universality of particle production at high energy. In addi- 
tion, the AuAu data approximately agrees with the scaled pp and pp data 
suggesting that the effective energy of a high energy AA collision is just 
~ ,ys NN . This result is not understood theoretically and remains a puzzle. 

6. The Latest Results from RHIC 

At the time of the lectures, the critical dAu "control" run at RHIC was not 
complete. Since then, results from this run have been published by all four 
collaborations 34 > 57 > 58 ' 59 . These results show no hadron or jet suppression 
in dAu implying that the suppression is NOT present in the nuclear wave- 
function. This strongly favors the jet quenching interpretation for hadron 
suppression in AuAu and has led to a lot of theoretical activity. 




Figure 23. The nuclear modification factor. Left panel: Midrapidity result from 
PHENIX. Minimum bias dAu charged hadron result compared to central AuAu charged 
hadrons and minimum bias dAu neutral pions. 57 . Right panel: PHOBOS results slightly 
forward of mid-rapidity (0.2 < r\ < 1.4). The centrality dependence of the nuclear mod- 
ification factor for charged hadrons in dAu compared to central AuAu 58 . 

Figure 23 (left panel) shows the minimum bias dAu results from 
PHENIX at midrapidity. The charged hadrons are enhanced rather than 
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suppressed, in sharp contrast to AuAu, a result confirmed by BRAHMS 34 . 
Furthermore, the pions show collision scaling, again in sharp contrast to 
the strong suppression seen in AuAu. 

This contrast is striking, but the comparison of minimum bias dAu to 
central AuAu is not fully decisive. Any nuclear effects in dAu are expected 
to manifest themselves primarily in central collisions and can be washed 
out in minimum bias collisions. Fig. 23 (right panel) shows the centrality 
dependence of RdAu from PHOBOS, slightly forward of midrapidity (from 
the deuteron's point of view). Even the most central dAu collisions show 
no suppression. 




n Aty (radians) 



Figure 24. Azimuthal correlations for minimum bias dAu, central dAu, and central 
AuAu collisions compared to pp from STAR 59 . 



Finally, Fig. 24, from STAR, shows that the jet structure in central dAu 
collisions can be understood based on a pp reference sample. There is no 
significant reduction of back-to-back jets in head-on dAu collisions. The 
complete suppression of the away-side jet in central AuAu collisions is also 
repeated in this plot for comparison. 

Taken together, these results indicate that jets are quenched in AuAu 
collisions at RHIC energies (s/s NN of 130 and 200 GeV), while there is little 
or no evidence of such quenching in dAu or in lower energy A A. 
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The production and behavior of protons in dAu collisions is again sur- 
prising, however. Figure 23 shows that midrapidity pions scale like N co u 
at high momentum while total charged particles (including protons) are 
enhanced. This may explain why PHOBOS (Fig. 23) sees little enhance- 
ment of charged particles (fewer protons for rj ~ 0.8). The mystery comes 
from the fact that the explanations put forward for the relative behavior of 
protons and pions in AuAu do not explain their behavior in dAu. 

7. Summary 

The held of heavy ion physics has indeed reached an exciting new era. We 
have created a high temperature, high density, strongly interacting bulk 
matter state in the laboratory, and we have achieved temperatures higher 
than needed to theoretically create a quark-gluon plasma. This bulk matter 
exhibits interesting properties. It appears to be very dense and opaque even 
at high p T , generating the maximum possible elliptic flow and strongly 
quenching any jets which are not formed on the surface of the material. 
Furthermore, the system appears to be exploding in all three dimensions. 

Some puzzles remain. Why are there so many protons at high p T , and 
why do protons and pions behave differently even in dAu collisions? Is the 
particle production universal between AA, pp, pp, and e + e~ at high energy, 
and if so, why? 

Much work remains to be done to study this strongly interacting matter 
more quantitatively and to resolve the puzzles. Fortunately, the detectors 
and accelerator are undergoing continuous upgrades and the prospects for 
a continued rich harvest of physics from RHIC look excellent. 
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